Naval  Research  Laboratory 

Washington,  DC  20375-5320 


NRL/MR/6750-13-9489 


Finite  Larmor  Radius  and 
Collisional  Effects  on  the 
Electron-Ion  Hybrid  Instability 


E.M.  Tejero 
G.  Ganguli 
C.  Crabtree 
W.E.  Amatucci 

Charged  Particle  Physics  Branch 
Plasma  Physics  Division 


V.  Sotnikov 

Air  Force  Research  Laboratory 
Dayton,  Ohio 


C.L.  Enloe 

Air  Force  Academy 
Colorado  Springs,  Colorado 


November  18,  2013 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 
Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 
18-11-2013 


2.  REPORT  TYPE 

Memorandum  Report 


3.  DATES  COVERED  (From  -  To) 


4.  TITLE  AND  SUBTITLE 


Finite  Larmor  Radius  and  Collisional  Effects  on  the  Electron-Ion  Hybrid  Instability 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

E.M.  Tejero,  V.  Sotnikov,1  G.  Ganguli,  C.  Crabtree,  W.E.  Amatucci,  and  C.L.  Enloe2 


5d.  PROJECT  NUMBER 

67-9872-03  &  67-4774-02 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 
4555  Overlook  Avenue,  SW 
Washington,  DC  20375-5320 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


NRL/MR/6750— 13-9489 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


Office  of  Naval  Research 

One  Liberty  Center 

875  North  Randolph  Street,  Suite  1425 

Arlington,  VA  22203-1995 


Air  Force  Office  of  Scientific  Research 

875  North  Randolph  Street,  Suite  325,  Rm.  3112 

Arlington,  VA  22203-1768 


10.  SPONSOR  /  MONITOR’S  ACRONYM(S) 

ONR  &  AFOSR 


11.  SPONSOR  /  MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

1  Air  Force  Research  Laboratory,  Dayton,  OH 
2Air  Force  Academy,  Colorado  Springs,  CO 


14.  ABSTRACT 

The  Electron-Ion  Hybrid  instability,  a  transverse  velocity  shear-driven  instability  with  frequency  near  the  lower  hybrid  frequency,  has  been 
observed  theoretically  and  experimentally.  It  was  shown  previously  that  the  scale  length  of  the  gradient  in  the  velocity  must  be  much  smaller  than 
the  ion  gyroradius  and  larger  than  the  electron  gyroradius  in  order  to  generate  the  short  wavelength  electron-ion  hybrid  mode.  In  this  paper,  the 
original  theory  for  the  electron-ion  hybrid  instability  has  been  extended  to  include  finite  gyroradius  radius  effects  and  electron-neutral  collisions 
with  the  intention  of  applying  this  theory  to  the  plasma  region  surrounding  hypersonic  vehicles.  In  this  plasma  layer,  these  sorts  of  transverse 
sheared  flows  can  exist  in  a  collisional  plasma.  While  this  dense  layer  of  plasma  can  itself  impede  communications,  the  density  structures  created 
by  the  lower  hybrid  turbulence  can  also  be  a  source  of  scattering  for  these  electromagnetic  signals. 
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I.  INTRODUCTION 


Over  the  past  twenty  years,  there  have  been  numerous  studies  both  theoretical1-9  and 
experimental10-15  into  the  destabilizing  effects  of  inhomogeneous  electric  fields  transverse  to 
the  ambient  magnetic  held.  It  has  been  demonstrated  that  the  electric  held  gradient  scale 
length  Le  is  the  key  length  scale  that  determines  the  nature  of  the  instability16.  If  Le  is  on 
the  order  of  the  ion  gyroradius  pz ,  then  the  instability  will  be  driven  by  the  dynamics  of  the 
ions  and  could  be  an  electrostatic  ion  cyclotron  wave  or  an  Alfven  wave  depending  on  the 
plasma  /315.  However,  if  pe  <C  Le  <  Pi,  then  the  ions  do  not  ExB  drift  and  the  instability 
will  be  driven  by  the  dynamics  of  the  electrons.  These  waves  are  primarily  electrostatic 
lower  hybrid  waves  in  the  VLF  range.  All  of  these  velocity  shear-driven  instabilities  are 
primarily  azimuthally  propagating  in  the  direction  of  the  ExB  drift. 


In  this  report,  we  will  focus  on  the  electron-ion  hybrid  instability  (EIH),  which  is  observed 
in  the  case  where  pe  Le  <  Pi-  This  case  has  been  studied  in  relation  to  boundary  layer 
conditions  in  the  plasma  sheet  and  in  laser  generated  plasmas,  where  lower  hybrid  waves 
have  been  observed  in  environments  with  sharp  electric  held  gradients  that  are  required  for 
the  instability.  We  want  to  extend  this  work  to  include  the  plasma  environment  created 
around  hypersonic  vehicles.  Depending  on  the  speed  of  the  vehicle,  the  plasma  densities  can 
become  quite  large  such  that  the  plasma  can  cutoff  any  communications  signals  to  and  from 
the  vehicle.  However,  we  are  interested  in  densities  below  the  point  where  communication 
blackout  occurs.  The  plasma  layer  that  is  created  also  has  a  velocity  profile  due  to  motion 
of  the  vehicle  through  the  neutral  medium.  This  sheared  plasma  flow  can  give  rise  to 
similar  instabilities  to  those  described  above.  These  plasma  instabilities  in  their  nonlinearly 
saturated  state  can  form  coherent  density  structures  that  can  also  disrupt  communication 
signals.  It  is  these  nonlinear  structures  and  their  effects  on  electromagnetic  radiation  that  we 
wish  to  study.  To  that  end,  we  will  extend  the  EIH  theory  to  include  finite  electron-neutral 
collisions  and  finite  larmor  radius  effects  to  include  important  facets  of  the  hypersonic  vehicle 
plasma  environment. 

Manuscript  approved  September  10,  2013. 
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II.  THEORY 


In  order  to  extend  the  existing  theory  for  the  EIH  instability,  we  start  the  derivation  from 
the  fluid  equations  for  a  species  a  with  no  temperature  gradient: 

d  na 


dt 


+  V  ■  (nava)  =  0, 


dv. 


Qa 


0,  +  (va  •  v)  vQ  = — (E  +  vaxB) 
dt  mn 


T 

1  rv 


mana 


-Vn„. 


(1) 


(2) 


We  linearize  the  equations  assuming  A  ~  A0  +  Ai(x)  exp  [i(kyy  +  kzz  —  cot)}  +  ...  higher  order 
terms  that  we  neglect.  We  take  the  inhomogeneity  to  be  in  the  x  direction:  E  =  E0(a;)a;+E1, 
the  magnetic  held  to  be  B  =  B0z,  and  the  equilibrium  drift  along  v0a  =  v0ay. 

From  the  continuity  equation  to  first  order,  we  arrive  at  the  following  expressions: 
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^  ~dx  kzv\c 


(3) 


where  we  have  assumed  that  the  equilibrium  electron  and  ion  densities  are  uniform  and 
equal  to  n0  and  defined  coia  =  to  —  kyv0a. 

From  the  momentum  equation  to  first  order,  we  arrive  at  the  following  expression: 
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where  we  have  defined  Da  =  ui2a-  r]aQ2a,  cda  =  ujla  +  ivan,  rja  =  1  +  Aa  =  and 


vta  =  ^ 


Substituting  Equations  (4-6)  into  Equation  (3),  we  can  write  a  general  expression  for  the 
perturbed  density  for  a  species  a. 
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We  assume  that  the  ions  are  an  unmagnetized  cold  fluid  and  that  ion-neutral  collisions 
can  be  ignored.  These  assumptions  imply  that  the  ions  do  not  ExB  drift  v0i  =  0,  rji  —  1, 


2 


Uu  =  ca,  vin  =  0,  and  Tt  =  0.  In  addition,  the  frequency  of  interest  ca  is  near  the  lower  hybrid 
frequency  such  that  fle  S>  ta  S>  fl*.  This  implies  that  Dx  «  c a2.  Under  these  assumptions, 
we  can  rewrite  Equation  (7)  for  the  ions  as: 


nu  = 


en0 


rrijOJ 2 


d'2 

_ _ k2  -  k2 

dx2  y  z 


(8) 


The  electrons  are  taken  to  be  a  warm  magnetized  fluid,  and  we  keep  the  electron-neutral 
collision  frequency.  The  electrons  ExB  drift  at  Voe  =  vE  =  —  j^-  As  for  the  ions  we  treat 
the  frequency  such  that  f2e  3>  ca  3>  fb,  which  implies  that  De  — 7/ef22.  We  also  consider 
these  equations  in  the  small  shear  limit  such  that  <C  1.  Under  these  assumptions,  we 

can  rewrite  Equation  (7)  for  the  electrons  as: 
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The  partial  x  derivative  has  a  scale  size  equal  to  the  electric  held  gradient  scale  length  LE- 
Since  pe  <C  Le  <  Pi,  which  is  why  the  electrons  ExB  drift  and  the  ions  do  not,  we  can  drop 
the  second  derivative  of  n \ e  since  it  goes  like  p2/ L2E.  This  leads  to  the  following  expression 
for  the  perturbed  electron  density: 
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where  k2  =  ^1  +  (k2  — 
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We  substitute  Equation  (8)  and  Equation  (10)  into  Poisson’s  equation,  V20i  +  Wq 


p^'Hie  =  0,  and  after  a  mess  of  algebra,  we  arrive  at: 
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where  <5  =  cape/fle,  A2  =  2(Up2/(l  +  (r),  M  =  rrii/rne  and  we  have  ignored  terms  with  factors 


of  venR2- 


III.  NUMERICAL  SOLUTIONS 

Equation  (11)  is  an  eigenvalue  equation  that  can  be  solved  numerically  for  the  wave  po¬ 
tential  profile  and  the  complex  eigenvalue  ca.  We  can  rewrite  Equation  (11)  using  normalized 


3 


quantities  as  follows: 


where  we  have  defined  Q  as 
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1  +  i'2 

and  we  have  used  the  lower  hybrid  frequency  c olh  and  the  electric  held  scale  length  LE  to 
dehne  our  characteristic  time  and  length  scales  for  the  normalized  quantities,  e.g.  ky  =  kyLE 
and  Cu  —  uj/ljlh-  It  follows  that  the  characteristic  velocity  scale  is  ujlhLe  such  that  vE  = 
ve/(wlhLe)-  It  should  also  be  noted  that 
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For  onr  numerical  calculations,  we  will  be  using  the  following  expression  for  R2  that  is  derived 
from  a  kinetic  treatment17  instead  of  the  expression  determined  in  the  previous  section  from 
the  fluid  approximation: 


for  !  L  3>  u. 


R2  = 


pe 


3 

-Pg  +  3/9  ■ ,  for  ne  <C  CUpe 


(15) 


where  the  ion  Debye  length  is  \2Di  =  R/ (Anne2). 

If  we  consider  large  x  values,  where  we  can  take  vE  and  its  derivatives  to  be  zero,  Equation 
(13)  has  no  spatial  dependence  and  Equation  (12)  can  be  immediately  solved  to  yield: 


(x)  =  A  exp  [kxx)  +  B  exp  (—kxx), 


(16) 


where  kx  =  ±a JQ  (x  — >  ±oo,ta).  To  start,  we  choose  the  same  electric  held  profile  used  i 
Ganguli,  Lee,  and  Palmadesso2: 

E0 


m 


Ex  = 


cosh2  (x) 


-x. 


(17) 


We  use  a  standard  numerical  shooting  code18  to  solve  Equation  (12).  Using  the  asymp¬ 
totic  solution  given  in  Equation  (16),  we  start  the  the  numerical  integration  with  a  purely 
decaying  solution  for  a  large  negative  value  of  x  for  a  given  guess  for  uj.  We  determine  the 
values  of  the  coefficients,  A  and  B ,  from  the  result  of  the  numerical  integration  for  a  large 
positive  value  of  x.  We  adjust  Co  to  minimize  the  coefficient  in  front  of  the  growing  asymp¬ 
totic  solution.  Figure  (1)  (a)  shows  the  real  (solid  black  line)  and  imaginary  (dashed  red 
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(a) 


(b) 


FIG.  1.  Plot  (a)  shows  the  real  (solid  black  line)  and  imaginary  (dashed  red  line)  parts  of  a  typical 
eigenfunction  solution  of  (12)  with  kyLE  =  1.0,  5  =  10.0,  a  =  0 .1,  and  u  =  (1.9  +  0.94 i)ujlh-  Plot 
(b)  shows  profiles  of  the  real  frequency  (solid  black  lines)  and  growth  rate  (dashed  red  lines)  as 
functions  of  kyLE  for  three  different  values  of  5. 


line)  parts  of  a  typical  eigenfunction  solution  of  (12)  for  the  following  parameters:  ky  =  1.0, 
kz  =  0,  5  =  10.0,  ve  =  0.1005,  M  =  1833,  uen  =  0,  and  pe  =  0.  The  resulting  complex 
eigenvalue  is  to  =  (1.9  +  0.9Ai)uiLH-  Figure  (1)  (b)  shows  profiles  of  the  real  frequency  (solid 
black  lines)  and  growth  rate  (dashed  red  lines)  as  functions  of  ky  for  three  different  values 
of  5.  These  results  are  consistent  with  figures  in  the  Ganguli,  Lee,  and  Palmadesso  paper2, 
which  were  determined  for  similar  parameters.  With  these  benchmarks  for  the  shooting 
code,  we  can  now  look  at  how  ven  and  pe  effect  these  solutions. 


IV.  DISCUSSION 

The  plasma  environment  near  hypersonic  vehicles  can  be  highly  collisional,  so  we  wish  to 
investigate  the  effects  of  neutral  collisions  on  the  electron-ion  hybrid  instability.  Since  the 
ions  are  already  considered  as  being  unmagnetized,  the  effects  of  ion-neutral  collisions  will 
be  minimal,  and  consequently  we  ignored  them  in  the  above  derivation.  The  frequency  of 
the  EIH  instability  tends  to  be  near  the  lower  hybrid  frequency  though  it  can  be  Doppler- 
shifted  significantly  due  to  the  presence  of  the  nonuniform  electron  flows.  We  expect  to 
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FIG.  2.  Plot  shows  the  growth  rate  as  a  function  of  normalized  electron-neutral  collision  frequency 
for  three  values  of  6  and  the  same  conditions  used  in  Figure  (lb). 


start  to  see  the  effects  of  the  electron-neutral  collisions  when  uen  «  ujlh-  Due  to  the  drift 
approximation  that  we  made  in  the  above  derivation,  there  is  an  upper  limit  to  the  collision 
frequency  we  can  study  with  Equation  (12),  ven  <  0.1  f2e  such  that  uen/Qe  <C  1.  Figure  (2) 
shows  typical  profiles  of  growth  rate  as  a  function  of  the  electron- neutral  collision  frequency 
normalized  to  the  lower  hybrid  frequency  for  the  same  three  values  of  5  used  in  Figure  (lb), 
where  the  range  of  the  abscissa  is  [10-5  fle,  0.1  fle].  As  can  be  seen  there  is  a  sudden 
drop  in  the  growth  rate  for  the  three  curves,  but  the  plasma  is  not  completely  stabilized 
by  electron-neutral  collisions.  The  critical  values  where  collisional  effects  become  important 
are  dependent  on  the  value  of  8:  for  5  =  10.0  =>-  ven  =  0.22  u>lh ,  5  =  1.0  =>■  uen  =  0.22  u>lh, 
and  5  =  0.5  =>■  uen  =  1.35  colh- 

We  seek  to  investigate  the  nonlinear  saturation  of  this  instability.  Specifically,  we  wish 
to  study  the  formation  of  coherent  stationary  structures  that  result  from  the  nonlinear 
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1.4 


FIG.  3.  Plot  shows  profiles  of  real  frequency  (solid  black  line)  and  growth  rate  (dashed  red  line) 
as  functions  of  kype.  These  profiles  were  calculated  using  the  following  parameters:  6  =  0.35, 
p  =  40.0,  a  =  0.002,  kz  =  uen  =  0,  and  pe  =  0.046. 


coupling  of  the  instability  in  the  lower  hybrid  range  and  ion  acoustic  waves  by  means  of 
the  modulation  instability.  In  order  to  have  this  nonlinear  pathway  be  dominant,  we  search 
for  solutions  that  are  centered  around  the  lower  hybrid  frequency  and  that  have  unstable 
solutions  for  a  narrow  range  of  kype.  Figure  (3)  shows  profiles  of  real  frequency  (solid  black 
line)  and  growth  rate  (dashed  red  line)  as  functions  of  kype.  These  profiles  were  calculated 
using  the  following  parameters:  S  =  0.35,  p  =  40.0,  a  =  0.002,  kz  =  ven  =  0,  and  pe  =  0.046. 
For  these  conditions,  the  unstable  solutions  are  centered  in  a  small  range  around  the  lower 
hybrid  frequency  and  occur  in  a  narrow  range  of  kype.  In  future  studies,  the  following 
nonlinear  equations  from  Sotnikov  et  ah19  will  be  solved  to  investigate  the  development  of 
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the  nonlinear  coherent  stationary  structures: 
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where  we  have  used  the  following  definitions 
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dxdy  dy  dx 


,  f  d  .  .  d 
and  L1  =  —+ve{x)—. 
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(18) 
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V.  CONCLUSION 

We  have  extended  the  work  of  Ganguli  et  al.  for  the  electron-ion  hybrid  in  order  to 
include  important  effects  for  the  plasma  environment  surrounding  a  hypersonic  vehicle  and 
written  the  necessary  shooting  code  to  solve  this  new  system.  We  have  shown  that  the 
instability  is  robust;  the  effects  of  electron- neutral  collisions  take  effect  when  uen  «  uJlh- 
However,  even  electron-neutral  collision  frequencies  near  0.1  Vte  do  not  fully  stabilize  the 
plasma,  and  above  this  threshold  the  notion  of  an  ExB  drift  is  no  longer  valid.  With 
finite  larmor  radius  effects  included,  we  have  obtained  an  operating  range  where  we  can  look 
for  the  development  of  nonlinear  stationary  structures  due  to  the  nonlinear  coupling  of  the 
unstable  shear-driven  lower  hybrid  mode  and  low  frequency  ion  acoustic  waves. 
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